Structure-based design and in-parallel synthesis of inhibitors of AmpC beta-lactamase

Background: Group I p-lactamases are a major cause of antibiotic resistance to beta -lactams such as penicillins and cephalosporins. These enzymes are only modestly affected by classic beta -lactam-based inhibitors, such as clavulanic acid. Conversely, small arylboronic acids inhibit these enzymes at sub-micromolar concentrations. Structural studies suggest these inhibitors bind to a well-defined cleft in the group I beta -lactamase AmpC; this cleft binds the ubiquitous R1 side chain of beta -lactams. Intriguingly, much of this cleft is left unoccupied by the small arylboronic acids. Results: To investigate if larger boronic acids might take advantage of this cleft, structure-guided in-parallel synthesis was used to explore new inhibitors of AmpC. Twenty-eight derivatives of the lead compound, 3-aminophenylboronic acid, led to an inhibitor with 80-fold better binding (2; K-i 83 nM). Molecular docking suggested orientations for this compound in the R1 cleft. Based on the docking results, 12 derivatives of 2 were synthesized, leading to inhibitors with iii values of 60 nM and with improved solubility. Several of these inhibitors reversed the resistance of nosocomial Gram-positive bacteria, though they showed little activity against Gram-negative bacteria. The X-ray crystal structure of compound 2 in complex with AmpC was subsequently determined to 2.1 Angstrom resolution. The placement of the proximal two-thirds of the inhibitor in the experimental structure corresponds with the docked structure, but a bond rotation leads to a distinctly different placement of the distal part of the inhibitor. In the experimental structure, the inhibitor interacts with conserved residues in the R1 cleft whose role in recognition has not been previously explored. Conclusions: Combining structure-based design with in-parallel synthesis allowed for the rapid exploration of inhibitor functionality in the R1 cleft of AmpC. The resulting inhibitors differ considerably from beta -lactams but nevertheless inhibit the enzyme well. The crystal structure of 2 (K-i 83 nM) in complex with AmpC may guide exploration of a highly conserved, largely unexplored cleft, providing a template for further design against AmpC beta -lactamase

Combinatorial Computational Approaches to Identify Tetracycline Derivatives as Flavivirus Inhibitors

Limited structural information of drug targets, cellular toxicity possessed by lead compounds, and large amounts of potential leads are the major issues facing the design-oriented approach of discovering new leads. In an attempt to tackle these issues, we have developed a process of virtual screening based on the observation that conformational rearrangements of the dengue virus envelope protein are essential for the mediation of viral entry into host cells via membrane fusion. Screening was based solely on the structural information of the Dengue virus envelope protein and was focused on a target site that is presumably important for the conformational rearrangements necessary for viral entry. To circumvent the issue of lead compound toxicity, we performed screening based on molecular docking using structural databases of medical compounds. To enhance the identification of hits, we further categorized and selected candidates according to their novel structural characteristics. Finally, the selected candidates were subjected to a biological validation assay to assess inhibition of Dengue virus propagation in mammalian host cells using a plaque formation assay. Among the 10 compounds examined, rolitetracycline and doxycycline significantly inhibited plaque formation, demonstrating their inhibitory effect on dengue virus propagation. Both compounds were tetracycline derivatives with IC(50)s estimated to be 67.1 µM and 55.6 µM, respectively. Their docked conformations displayed common hydrophobic interactions with critical residues that affected membrane fusion during viral entry. These interactions will therefore position the tetracyclic ring moieties of both inhibitors to bind firmly to the target and, subsequently, disrupt conformational rearrangement and block viral entry. This process can be applied to other drug targets in which conformational rearrangement is critical to function

The first small-molecule inhibitors of members of the ribonuclease E family

The Escherichia coli endoribonuclease RNase E is central to the processing and degradation of all types of RNA and as such is a pleotropic regulator of gene expression. It is essential for growth and was one of the first examples of an endonuclease that can recognise the 5′-monophosphorylated ends of RNA thereby increasing the efficiency of many cleavages. Homologues of RNase E can be found in many bacterial families including important pathogens, but no homologues have been identified in humans or animals. RNase E represents a potential target for the development of new antibiotics to combat the growing number of bacteria that are resistant to antibiotics in use currently. Potent small molecule inhibitors that bind the active site of essential enzymes are proving to be a source of potential drug leads and tools to dissect function through chemical genetics. Here we report the use of virtual high-throughput screening to obtain small molecules predicted to bind at sites in the N-terminal catalytic half of RNase E. We show that these compounds are able to bind with specificity and inhibit catalysis of Escherichia coli and Mycobacterium tuberculosis RNase E and also inhibit the activity of RNase G, a paralogue of RNase E

Core Site-Moiety Maps Reveal Inhibitors and Binding Mechanisms of Orthologous Proteins by Screening Compound Libraries

Members of protein families often share conserved structural subsites for interaction with chemically similar moieties despite low sequence identity. We propose a core site-moiety map of multiple proteins (called CoreSiMMap) to discover inhibitors and mechanisms by profiling subsite-moiety interactions of immense screening compounds. The consensus anchor, the subsite-moiety interactions with statistical significance, of a CoreSiMMap can be regarded as a “hot spot” that represents the conserved binding environments involved in biological functions. Here, we derive the CoreSiMMap with six consensus anchors and identify six inhibitors (IC50<8.0 µM) of shikimate kinases (SKs) of Mycobacterium tuberculosis and Helicobacter pylori from the NCI database (236,962 compounds). Studies of site-directed mutagenesis and analogues reveal that these conserved interacting residues and moieties contribute to pocket-moiety interaction spots and biological functions. These results reveal that our multi-target screening strategy and the CoreSiMMap can increase the accuracy of screening in the identification of novel inhibitors and subsite-moiety environments for elucidating the binding mechanisms of targets

Discovery of novel GPVI receptor antagonists by structure-based repurposing.

Inappropriate platelet aggregation creates a cardiovascular risk that is largely managed with thienopyridines and aspirin. Although effective, these drugs carry risks of increased bleeding and drug 'resistance', underpinning a drive for new antiplatelet agents. To discover such drugs, one strategy is to identify a suitable druggable target and then find small molecules that modulate it. A good and unexploited target is the platelet collagen receptor, GPVI, which promotes thrombus formation. To identify inhibitors of GPVI that are safe and bioavailable, we docked a FDA-approved drug library into the GPVI collagen-binding site in silico. We now report that losartan and cinanserin inhibit GPVI-mediated platelet activation in a selective, competitive and dose-dependent manner. This mechanism of action likely underpins the cardioprotective effects of losartan that could not be ascribed to its antihypertensive effects. We have, therefore, identified small molecule inhibitors of GPVI-mediated platelet activation, and also demonstrated the utility of structure-based repurposing

Multilevel Parallelization of AutoDock 4.2

<p>Abstract</p> <p>Background</p> <p>Virtual (computational) screening is an increasingly important tool for drug discovery. AutoDock is a popular open-source application for performing molecular docking, the prediction of ligand-receptor interactions. AutoDock is a serial application, though several previous efforts have parallelized various aspects of the program. In this paper, we report on a multi-level parallelization of AutoDock 4.2 (mpAD4).</p> <p>Results</p> <p>Using MPI and OpenMP, AutoDock 4.2 was parallelized for use on MPI-enabled systems and to multithread the execution of individual docking jobs. In addition, code was implemented to reduce input/output (I/O) traffic by reusing grid maps at each node from docking to docking. Performance of mpAD4 was examined on two multiprocessor computers.</p> <p>Conclusions</p> <p>Using MPI with OpenMP multithreading, mpAD4 scales with near linearity on the multiprocessor systems tested. In situations where I/O is limiting, reuse of grid maps reduces both system I/O and overall screening time. Multithreading of AutoDock's Lamarkian Genetic Algorithm with OpenMP increases the speed of execution of individual docking jobs, and when combined with MPI parallelization can significantly reduce the execution time of virtual screens. This work is significant in that mpAD4 speeds the execution of certain molecular docking workloads and allows the user to optimize the degree of system-level (MPI) and node-level (OpenMP) parallelization to best fit both workloads and computational resources.</p

Specificity quantification of biomolecular recognition and its implication for drug discovery

Highly efficient and specific biomolecular recognition requires both affinity and specificity. Previous quantitative descriptions of biomolecular recognition were mostly driven by improving the affinity prediction, but lack of quantification of specificity. We developed a novel method SPA (SPecificity and Affinity) based on our funneled energy landscape theory. The strategy is to simultaneously optimize the quantified specificity of the “native” protein-ligand complex discriminating against “non-native” binding modes and the affinity prediction. The benchmark testing of SPA shows the best performance against 16 other popular scoring functions in industry and academia on both prediction of binding affinity and “native” binding pose. For the target COX-2 of nonsteroidal anti-inflammatory drugs, SPA successfully discriminates the drugs from the diversity set, and the selective drugs from non-selective drugs. The remarkable performance demonstrates that SPA has significant potential applications in identifying lead compounds for drug discovery

Intramolecular Epistasis and the Evolution of a New Enzymatic Function

Atrazine chlorohydrolase (AtzA) and its close relative melamine deaminase (TriA) differ by just nine amino acid substitutions but have distinct catalytic activities. Together, they offer an informative model system to study the molecular processes that underpin the emergence of new enzymatic function. Here we have constructed the potential evolutionary trajectories between AtzA and TriA, and characterized the catalytic activities and biophysical properties of the intermediates along those trajectories. The order in which the nine amino acid substitutions that separate the enzymes could be introduced to either enzyme, while maintaining significant catalytic activity, was dictated by epistatic interactions, principally between three amino acids within the active site: namely, S331C, N328D and F84L. The mechanistic basis for the epistatic relationships is consistent with a model for the catalytic mechanisms in which protonation is required for hydrolysis of melamine, but not atrazine

A statistical framework to evaluate virtual screening

<p>Abstract</p> <p>Background</p> <p>Receiver operating characteristic (ROC) curve is widely used to evaluate virtual screening (VS) studies. However, the method fails to address the "early recognition" problem specific to VS. Although many other metrics, such as RIE, BEDROC, and pROC that emphasize "early recognition" have been proposed, there are no rigorous statistical guidelines for determining the thresholds and performing significance tests. Also no comparisons have been made between these metrics under a statistical framework to better understand their performances.</p> <p>Results</p> <p>We have proposed a statistical framework to evaluate VS studies by which the threshold to determine whether a ranking method is better than random ranking can be derived by bootstrap simulations and 2 ranking methods can be compared by permutation test. We found that different metrics emphasize "early recognition" differently. BEDROC and RIE are 2 statistically equivalent metrics. Our newly proposed metric SLR is superior to pROC. Through extensive simulations, we observed a "seesaw effect" – overemphasizing early recognition reduces the statistical power of a metric to detect true early recognitions.</p> <p>Conclusion</p> <p>The statistical framework developed and tested by us is applicable to any other metric as well, even if their exact distribution is unknown. Under this framework, a threshold can be easily selected according to a pre-specified type I error rate and statistical comparisons between 2 ranking methods becomes possible. The theoretical null distribution of SLR metric is available so that the threshold of SLR can be exactly determined without resorting to bootstrap simulations, which makes it easy to use in practical virtual screening studies.</p